The depletion force in a bi-disperse granular layer 
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We demonstrate the effect of the depletion force in experiments and simulations of vertically 
vibrated mixtures of large and small steel spheres. The system exhibits size segregation and a 
large increase in the pair correlation function of the large spheres for short distances that can be 
accurately described using a combination of the depletion potential derived for equilibrium colloidal 
systems and a Boltzmann factor. The Boltzmann factor defines an effective temperature for the 
system, which we compare to other measures of the temperature. 
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Vibrated granular materials show many complex be- 
haviors, including pattern formation, segregation, non- 
equipartition of energy, and non- Gaussian velocity dis- 
tributions Q, 0, • Many of these are purely nonequilib- 
rium effects, but some features of vibrated granular ma- 
terials are strikingly similar to equilibrium systems such 
as colloidal suspensions 0, 0, 0. The extent to which 
the machinery of equilibrium statistical mechanics can 
be used as a starting point to develop a theory of these 
non-equilibrium steady states remains an open question. 

Segregation mechanisms, meanwhile, are important in 
many industrial applications. Granular mixtures are 
known to segregate by size, mass, shape, and frictional 
coefficients and this segregation can be caused by many 
mechanisms, including vibration, convection, and tum- 
bling 0, 0- In this Letter, we investigate a relatively 
unexplored mechanism for segregation in granular mate- 
rials: the depletion force. Well studied in colloidal sys- 
tems H El E3 , the depletion force is an entropic effect 
in which the introduction of small particles into a col- 
loidal dispersion will cause the larger colloidal particles 
to segregate. Previous suggestions of the role of the de- 
pletion force in granular materials have focused on hor- 
izontally shaken [ill . HiL or horizontally swirled [l^ 
mixtures. The segregation effects for the horizontally vi- 
brated mixtures, however, has been shown to be caused 
by a differential driving mechanism which may also 
hold for the horizontally swirled system. In this paper 
we show that the depletion force can cause size segrega- 
tion in a vertically vibrated, bi-disperse granular layer. 
The pair correlation function for low concentrations of 
the large particles can be accurately described using the 
depletion potential derived for equilibrium systems and 
a Boltzmann factor. This Boltzmann factor defines an 
effective temperature for the system, which we compare 
to other measures of temperature in the system. 

The depletion force was originally proposed by 
Asakura and Oosawa in 1954 There are two ap- 
proaches for deriving the potential: the first is thermo- 
dynamic and the second is mechanical. In the thermody- 
namic view, illustrated in Fig. d the depletion force is 
derived through entropy maximization. There is a shell of 
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FIG. 1: Schematic drawing of the depletion interaction. The 
white shell around the large particles represents the volume 
from which the small particles are excluded. The dark blue 
region shows the overlap in this excluded volume when the 
large particles are close together. 



width r s (the radius of the smaller particles) around the 
large spheres into which the centers of the small spheres 
cannot penetrate. By arranging two large spheres to- 
gether, the excluded volume regions of the two large 
spheres overlap, and the volume accessible to the small 
spheres increases, increasing the total entropy of the sys- 
tem. This produces an effective attraction between the 
large spheres. The mechanical view of the depletion force 
is that the small spheres exert a pressure on the large 
spheres, but that pressure is unbalanced when two large 
spheres are close enough together that the small spheres 
are excluded from the space between the spheres. The 
net force is found by integrating the pressure over the 
surface of the large spheres. In equilibrium, both views 
are equivalent and give the following expression for the 
depletion potential for two large spheres separated by a 
distance r, valid for moderately low densities of small 
spheres: 
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where vl is the radius of the large spheres, P is the pres- 
sure exerted on the large spheres by the small spheres 
and V ex is the overlap in excluded volume. V ex can be 
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expressed as: 
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Our experimental setup consists of steel spheres of two 
sizes placed between two horizontal plates which are vi- 
brated sinusoidally in the vertical direction. Images cap- 
tured by a high resolution digital camera are analyzed 
to find the positions of the large spheres. For a range of 
ball diameters, as well as a range of both large and small 
ball volume fractions, we observe segregation of the large 
spheres. Fig. [2t shows the results of one such experi- 
ment, where large spheres initially distributed randomly 
on the plate have formed a compact cluster. The pa- 
rameter values, including the gap between the bottom 
and top plates, ft, the vibration frequency, the dimen- 
sionless acceleration T = A(2ttu) 2 jg, where A is the vi- 
bration amplitude, and the volume fraction of large and 
small spheres, pL, s — Nl, s Vl, s /V \ where Vl, s is the vol- 
ume of the large or small spheres and V is the volume of 
the system, are given in the figure caption. The segrega- 
tion behavior is relatively insensitive to the frequency and 
amplitude of the acceleration, but is sensitive to the gap 
spacing. For spacings very close to the large ball diame- 
ter, friction with the top and bottom plates severely lim- 
its motion of the large spheres. For spacings significantly 
larger than the large sphere diameter, the kinetic energy 
of the large spheres is much larger than that of the small 
spheres, and the forces from the small spheres do not ap- 
preciably affect the trajectories of the large spheres. This 
strong deviation from equipartition is a consequence of 
the strong interaction with the vibrating plate experi- 
enced by all of the particles in the thin granular layers. 
Except as noted, all of the results presented below were 
taken with the parameters used for Fig. [3 

Because the system is highly sensitive to variations in 
the flatness of the plates, we also performed molecular 
dynamics simulations of the system to ensure that the 
segregation was not caused by such experimental issues, 
using a model that has been shown to reproduce the dy- 
namics of similar granular experiments \% [ly, . The 
coefficient of restitution used in the simulation, which is 
set by a combination of the elastic restoring force and 
the dissipation, was 0.9. The other parameters were cho- 
sen to match the experiment. The simulations produced 
phase segregation similar to that observed in the experi- 
ment, one example of which is shown in Fig. 

In equilibrium hard sphere systems, the depletion force 
can be measured directly by looking at the pair correla- 
tion function. If the density of large spheres is small 
enough that only two-particle interactions are present, 
then the pair correlation function, G(r), can be sim- 
ply described through the use of a Boltzmann factor: 
G(r) = exp(— U(r)/T), in units where ks = 1- In 
the non-equilibrium granular system, there is no reason 




FIG. 2: Segregation due to the depletion force. Left panel: 
Experiment (p s = 0.036, pL = 0.066). Right Panel: Sim- 
ulations (p s = 0.080, pL = 0.132). In both instances the 
large and small spheres were initially randomly distributed 
throughout the available space. z/ = 30 Hz, r = 3.5, tl = 4 
mm, rz,/r s =3.33, ft — 2.3t\l. 



to expect that this equilibrium approach should accu- 
rately describe the pair correlation function. Fig. [Ht 
shows an example of the pair correlation function mea- 
sured in the experiment, and a result from the simula- 
tion is shown in Fig. I^Jd. The solid curve is a fit to 
G(r) = exp(-V ex P/T), with V ex given by Eq. Eland 
P/T a fitting parameter that provides only an overall 
scale factor. The equilibrium-like approach accurately 
reproduces the width of the peak in G(r) and the overall 
shape of the curve. We performed simulations with the 
small sphere radius reduced by 60% and found similar 
agreement. The fact that the range of the interaction 
between the large spheres is set by the small sphere ra- 
dius is unequivocal evidence that the attraction is due to 
the depletion force. 

In the simulations, the pressure in the system can 
be calculated by use of the virial theorem 0, EH, al- 
lowing an explicit calculation of an 'effective tempera- 
ture', Trfep, from the value of P/T obtained from the fit 
to G(r). In order to determine if this temperature is 
related to other fluctuations in the bi-disperse system, 
we have measured three other effective temperatures in 
the computer simulations, the granular temperature of 
the small spheres, the granular temperature of the large 
spheres, and an effective temperature defined through a 
fluctuation-dissipation relation. The granular tempera- 
tures are given by m < vi > 2 , where V{ is one component 
of the particle velocity and m is its mass. The distribu- 
tions of the horizontal velocity components for the large 
and small particles are shown in Fig. Both are ap- 
proximately Gaussian for small values of the velocity, and 
the small particles are well described by a gaussian dis- 
tribution over all of the accessible range. The smaller 
particles are less confined and hence have a more three- 
dimensional motion, which has been shown to lead to 
more nearly Gaussian velocity distributions [l8| . 

The granular temperatures are compared to Td ep for 
several concentrations of small spheres in Fig. 0t. The 
results show that, unlike either granular temperature, 
Td ep is a strong function of small sphere density. One 
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FIG. 3: Pair correlation function of large particles and the fit 
using the a Boltzmann factor and the equilibrium depletion 
potential (Eq. a) Experiment (p a = 0.071, p L = 0.007) b) 
Simulations (p s = 0.08, pL — 0.03) c) Logarithm of G(r) for 
several simulation runs (p s = 0.06, pL — 0.026). Solid lines 
are fits to the equilibrium theory (see text). 



possible explanation for this behavior is that the pres- 
sure in the vicinity of the large balls differs significantly 
from the average pressure due to non-equilibrium effects. 
We also investigated whether T& ev is related to an effec- 
tive temperature determined by a fluctuation-dissipation 
relation. While the granular system is not in the lin- 
ear non-equilibrium regime, some studies have found that 
fluctuation-dissipation measurements of the temperature 
can agree with the granular temperature [2(3, l-llj • In 
particular, in Ref. |2l[, it was found that the granular 
temperature of each species of particles in a binary gran- 
ular gas matches the temperature calculated by apply- 
ing the Einstein relation for that species. The Einstein 
relation connects the coefficient of viscous drag, &drag 5 
to the diffusion constant, D, through the temperature: 
T = kdragD. In the simulations, we can apply a con- 
stant external force to a large sphere and measure &drag 
using the relation F = ^drag^term, wnere ^term is the ter- 
minal velocity. An example of the linear displacement 
vs. time in the presence of a constant external force is 
shown in Fig. 2Jd. We can also measure the diffusion 
constant in the case of no external forces using the re- 



FIG. 4: Temperatures measured in the simulations of the 
granular system, a) Velocity distributions of small and large 
particles in the simulations with the same parameters as in 
Fig. 0: (vq = V< v 2 >). The solid line is a Gaussian dis- 
tribution, (b) Average displacement < x > as a function 
of time for a large particle subjected to a constant external 
force, and mean squared displacement < x 2 > vs. time for an 
unforced particle (p s = 0.096, pL = 0.0026). (c) The granu- 
lar temperatures of the small (black squares) and large (red 
circles) spheres and the fluctuation-dissipation temperature 
(green triangles) at different densities of the small spheres 
plotted vs. the effective temperature measured from for the 
depletion potential (From lowest to highest Td ep ' p s =0.036, 
0.060, 0.072, 0.097.) 

lation < x 2 >= 2Dt. An example of this is also shown 
in Fig. Et>. Combining these measurements gives us an 
independent measurement of temperature in the system: 
Tfd = DF ext /v term . As can be seen in Fig. Efc, the 
fluctuation-dissipation measurement of the temperature 
does not match the depletion temperature. Also, unlike 
the results described in Ref. [2l[, there is no agreement 
between the fluctuation-dissipation temperature and ei- 
ther granular temperature, although it is relatively close 
to the value found for the large spheres. This difference 
may be due to the forcing or dissipation from the top and 
bottom boundaries, which were not present in the model 
of Ref. [2l[. Similarly, in quasi-2D colloidal systems, 
the presence of the confining plates results in some small 
additional contributions to the depletion interaction 
that have not been considered in the above calculations. 
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As the density of small balls increases, the simple form 
of the depletion potential shown in Eq. is no longer 
accurate in equilibrium systems [2^. From an entropic 
point of view, the most favorable separations are those 
that allow integer numbers of small spheres to fit between 
large spheres. This leads to additional peaks in the pair 
correlation function in the equilibrium case at integer 
multiples of the small sphere diameter. This many-body 
effect, which is not easy to explain from the mechanical 
arguments, shows up as a small peak in G(r) for r = 
2vl + 2r s . This peak is visible in the pair correlation 
function determined from a simulation at a high density 
of small spheres, as shown in Fig. This demonstrates 
that subtle entropy- driven dynamics are observable in 
this far from equilibrium granular system. 
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FIG. 5: Pair correlation function for large spheres in the case 
of a high concentration of small spheres (p s = 0.12, pL = 
0.026). The smaller second peak in G(r) occurs at r = 2tl + 
2r s (indicated by the arrow), where exactly one small sphere 
fits in between the large spheres. 

In summary, we have found direct evidence for the de- 
pletion force in a vibrated granular system. We observe 
segregation of large particles induced by the presence 
of small particles, and an increase in the pair correla- 
tion function of the large particles at low concentrations 
which is consistent with the equilibrium depletion poten- 
tial. We also observe, for high concentrations of small 
particles, a second peak in the large particle pair cor- 
relation function that is characteristic of the depletion 
potential. Using an approach from equilibrium statisti- 
cal mechanics, we measure an effective temperature for 
the system using the depletion potential and a Boltz- 
mann factor. This effective temperature does not agree 
with the granular temperatures of either sized particle 
or with a measurement of the temperature with a simple 
fluctuation-dissipation measurement. These observations 
suggest that the driving mechanisms of the depletion in- 
teraction can operate in far from equilibrium granular 
systems, but that the equilibrium tools for quantifying 



the effects of the depletion force need modification to ac- 
count for nonequilibrium effects such as energy injection 
and dissipation. 
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